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The NASA Constellation program is utilizing Computational Fluid Dynamics (CFD) 
predictions for generating aerodynamic databases and design loads for the Ares I, Ares I-X, 
and Ares V launch vehicles and for aerodynamic databases for the Orion crew exploration 
vehicle and its launch abort system configuration. This effort presents several challenges to 
applied aerodynamicists due to complex geometries and flow physics, as well as from the 
juxtaposition of short schedule program requirements with high fidelity CFD simulations. 
NASA TetrUSS codes (GridTool/VGRID/USM3D) have been making extensive 
contributions in this effort. This paper will provide an overview of several enhancements 
made to the various elements of TetrUSS suite of codes. Representative TetrUSS solutions 
for selected Constellation program elements will be shown. Best practices guidelines and 
scripting developed for generating TetrUSS solutions in a production environment will also 
be described. 


I. Introduction 

NASA’s Constellation program^ was initiated to address the needs of U.S. aeeess to spaee in the post Spaee 
Shuttle era with the mission to the moon, Mars and beyond. Unlike the Spaee Shuttle, where both erew and eargo 
are launehed simultaneously, the proposed arehiteeture for the Constellation program ineludes two separate launeh 
vehieles: the Ares I for the erew and the Ares V for the eargo. The Ares I erew launeh vehiele (CLV) is a two-stage 
roeket that ineorporates a modified five-segment solid roeket booster from the Shuttle program while the upper stage 
is powered by the J-2X engine from the Apollo program. The major eomponents of the vehiele as well as the 
proposed nominal flight trajeetory are shown eoneeptually^ in figure 1. The major eomponents of Ares-I CLV 
eonsist of the Upper Stage (US) and the First Stage (FS) Solid Roeket Booster (SRB) that ineludes the frustum and 
the aft- skirt. Ares-I is the launeh vehiele for the Orion Crew Exploration Vehiele (CEV) and its Launeh Abort 
System (LAS) as also illustrated in figure 1 . 

The aerodynamie databases used for developing the Ares launeh vehieles and the Orion CEV are being produeed 
from a eombination of wind tunnel testing and eomputational fluid dynamies (CFD) simulations. Sinee these vehiele 
geometries have sueh simple shapes, the original assumption was that CFD eould serve as the primary souree of 
aerodynamie data with wind tunnel testing in an anehoring role. However, the initial CFD analysis on these simple 
shapes quiekly exposed eomplex flow fields primarily in assoeiation with blunt bases, protuberanee effeets, and 
elose proximity jet interaetions. The eomplex flow fields were also found to be ehallenging to model with CFD. 
Henee, the historieal role of the wind tunnel as the primary souree of aerodynamie data was restored, despite its own 
set of limitations and uneertainties. Even with CFD in a supporting role, thousands of solutions have been generated 
yielding an impressive amount of eomputational data for the Ares and Orion projeets. 
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Within the subsonic/transonic/supersonic speed range, the eomputational aerodynamie database for Constellation 
was primarily generated from the NASA TetrUSSVuSMBD^’^ and OVERFLOW^ flow solvers. Sinee eode-to-eode 
uneertainty was a key eoneern, some additional support was also provided from the NASA FUN3D^ eode. The 
NASA Carts D inviseid eode^ was also used heavily for sereening broad ranges of parameter spaee within the Orion 
projeet. For the Ares flight vehieles, the USM3D eode served as the primary eomputational data souree with 
supporting role from the OVERFFOW eode. For the Orion CEV and FAS, the OVERFFOW eode provided the 
majority of eomputational data with USM3D serving in a support role for quantifying eode-to-eode uneertainty. In 
both eases, the CFD teams worked seamlessly aeross the NASA eenters to deliver the needed data on time, while 
eollaborating often to solve the unforeseen ehallenges that surfaeed. 

Historieally, TetrUSS has been a researeh tool. However, Ares and Orion were elearly development projeets that 
demanded produetion-level effieieneies for data generation. Henee, many ehallenges were eneountered while 
quiekly transforming a researeh-oriented tool into a produetion-oriented tool in the midst of a fast-paeed database 
generation eampaign. Some of the ehallenges that were eneountered during the applieation of TetrUSS tools in the 
Ares and Orion projeets were: 

Early ehallenges: 

• Required grid sizes were mueh larger than USM3D eould aeeommodate in memory, 

• Solver ineffieieneies for large grids prohibited produetion-level throughput with USM3D, 

• How to quantify solution uneertainty, 

• Running and managing a large number of eases, 

• Post-proeessing and managing a massive amount of data. 

Computational ehallenges for Ares Projeet launeh vehieles: 

• Need for 10-15 solutions per day on grids with 50 to 100+ million tetrahedral eells, 

• Effieient generation of tetrahedral grids suitable for flight Reynolds number on slender launeh vehiele with 
large disparity of geometrie seale from small protuberanees, 

• Numerieal stiffness assoeiated with modeling flows at flight Reynolds number, 

• Extraeting line loads for struetural design, 

• Extraeting and organizing aerodynamie loads on numerous individual protuberanees for struetural design. 
Computational ehallenges for Orion Projeet FAS: 

• Farge uneertainty in axial foree due to strong sensitivities to grid resolution and turbulenee models, 

• Inadequate turbulenee models and gridding strategies for predieting strong jet interaetion effeets from abort 
motors, 

• Need for autonomous solution-adaptive grid eapability for resolving eomplex jet-plume flows, 

• Throughput requests on order of 400 FAS jet interaetion solutions per month. 

This paper presents an overview of how most of these ehallenges were sueeessfully resolved. The paper is 
organized with Seetion II presenting a broad overview of the TetrUSS system. Seetion III ineludes a diseussion of 
key enhaneements to the grid generation eapability that enabled the rapid generation of high-quality tetrahedral 
grids. A deseription of the solution-adaptive grid approaeh is also ineluded. Seetion IV deseribes enhaneements to 
the flow solver that overeame show stoppers for large applieations sueh as memory bottleneeks, slow run times, 
erippling robustness issues for high-Maeh number flows and startup of jet flows, inadequate jet boundary eondition 
eapability, and defieient turbulenee models for strong jet/aero- interaetion flows. The seetion V deseribes the 
proeesses devised for running and post-proeessing large numbers of eases. Finally, Seetion VI highlights some 
representative applieations for the Ares I, IX, and V launeh vehieles and the Orion Fauneh Abort System (FAS). 

II. Overview of TetrUSS Tool Set 

The NASA Tetrahedral Unstruetured Software System (TetrUSS) was developed during 1990s to provide a rapid 
aerodynamie analysis and design eapability to applied aerodynamieists. The system is eomprised of loosely 
integrated, user-friendly software that enables the applieation of advaneed Euler and Navier- Stokes tetrahedral finite 
volume teehnology to eomplex aerodynamie problems. A sehematie view of the TetrUSS system is presented in 
figure 2. The system eonsists of eomponent software for setting up geometrie surfaee definitions (GridTool), 
generating tetrahedral grids (VGRID), eomputing Euler and Navier-Stokes flow solutions (USM3D), and extraeting 
meaningful information from analysis of results (SimpleView). The system also allows for imposing design or 
aeroelastie shape ehanges by interfaeing with other eodes, sueh as CDISC^. An overview of the eapabilities of 
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TetrUSS system in early 2000 is presented in referenee 3. The salient features of the VGRID grid generator, 
USM3D flow solver, and various visualization as well as data analysis tools and utilities are presented below. 

A. Grid Generation 

VGRID is a tetrahedral grid generator based on the Advaneing Front Method^^ (AFM) for generation of the 
‘inviseid’ field eells and the Advaneing Layers Method^^ (ALM) for generation of thin-layered ‘viseous’ eells. Both 
teehniques are based on a marehing proeess in whieh tetrahedral eells grow on an initial front (triangular surfaee 
mesh) and gradually aeeumulate in the field around the subjeet geometry. Unlike the eonventional AFM, whieh 
introduees eells in the field around in a totally unstruetured manner, the ALM generates layers of thin tetrahedral 
eells in a more orderly fashion (one layer at a time) while maintaining the flexibility of AFM. Onee the advaneing 
front proeess is eompleted in VGRID, an additional post-proeessing step is required using POSTGRID to elose any 
open poekets and to improve grid quality. VGRID input files are generated by an interaetive geometry manipulation 
program, GridTool^^. It ean import surfaee definitions from IGES files (NURBS surfaees and NURBS eurves) and 
from PLOT3D files (point definitions) and manipulate them to define neeessary geometrie {surface patches) and 
grid-spaeing {sources) parameters. It uses OpenGL for 3D graphies. The graphieal interfaee is based on the Fast 
Light Toolkit. GridTool is available for Mae and Linux systems. 

B. Flow Solver 

The USM3D"^’^ eode is a parallelized tetrahedral eell-eentered, finite volume Navier-Stokes flow solver. The term 
eell eentered means that the finite volume flow equations are solved at the eentroid of eaeh tetrahedral eell. Inviseid 
flux quantities are eomputed aeross eaeh tetrahedral eell faee using various upwind sehemes. Spatial diseretization is 
aeeomplished by a novel reeonstruetion proeess, based on an analytieal formulation for eomputing solution gradients 
within tetrahedral eells. The solution is advaneed in time by a seeond-order Newton time step seheme^^, or to a 
steady-state eondition by an implieit baekward-Euler seheme. Several turbulenee models are available: the Spalart- 
Allmaras (SA) one-equation model^, the two-equation k-& turbulenee model, the Menter Shear Stress Transport 
(SST) two-equation model, and the nonlinear Algebraie Reynolds Stress Models (ARSM) of Girimaji and 
Shih/Zhu/Lumley^^. Detaehed Eddy Simulation (DES) has been implemented in all of the turbulenee models. A 
eapability to trip the flow at speeified loeations on aerodynamie surfaees has been implemented for the k-£, 
turbulenee model, but fully turbulent flow is assumed for the results to follow. USM3D has eapabilities for dynamie 
grid motion and overset grids 

C. Post Processing 

TetrUSS system has an interaetive menu-driven post-proeessing program. Simple View, whieh allows users to 
visualize unstruetured grid and flow solution details. SimpleView uses OpenGL for 3D graphies and is available for 
Mae and Linux systems. The Mae version uses the native Coeoa user interfaee, while the Linux version uses a 
Motif-based interfaee. It ean read grid files at any stage of the grid generation proeess and present either the entire 
grid generated or the surfaee mesh only. VGRID front files ean also be read, either in eonjunetion with the grid or in 
a stand-alone mode. Several diagnostie tools are available inside SimpleView, partieularly for the viseous layers of 
the grid. Users ean verify the presenee of viseous layers above any faee on a viseous surfaee and obtain a eount of 
the number of eells in the layers. Simpleview reads USM3D nodal solution file to display five eonserved flow 
variables, as well as a number of derived variables, sueh as pressure eoeffieient or entropy. It ean perform surfaee 
(line) and volume (planar) euts. Surfaee data and planar euts may be displayed as solid-shaded or eontour lines. It 
ean display oil flow and veloeity veetors on the aerodynamie surfaees. It ean also display veloeity veetors and 
streamlines in the flow-field and find iso-surfaees. 

Standalone eodes have also been developed to transform TetrUSS grid and solution files for visualization and 
post-proeessing using some of the major eommereial software, sueh as TeePlot® , FieldView®, and EnSight®. Other 
stand-alone post-proeessing tools have also been developed for the data analysis. For example, a utility to extraet 
seetional load distributions is available that will be deseribed later in more detail. A TetrUSS grid and solution 
eutter^^ utility for generalized sub-domain extraetion of over 50 user-speeified flow funetions has reeently been 
developed for effieient post proeessing of solutions. 

D. Utilities 

Several useful utility eodes are maintained through a eommon interfaee ealled ‘usgutiF. They are: 

1 . Mirroring of grid or VGRID input file, 

2. Extraeting eross-seetional geometry and flow field data, 
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3. Estimating grid spacing and stretching parameters used in VGRID for viscous grids, and 

4. Scaling, translating, rotating of grid and swapping of coordinates. 

Additional standalone utilities have also been developed to enhanee TetrUSS produetivity. For example, a utility 
is developed to interpolate solution from eoarse grid to a finer adapted grid. Another utility is available to identify 
trouble spots in difficult-to-eonverge or diverging solutions that ean aid in the localized improvements to the grid 
quality. 


III. Enhancements to Grid Generation 

Grid generation presents a number of challenges within a produetion environment. Geometries are continually 
evolving and require frequent generation of new grids. The demands of airframe design neeessitate adequate grid 
resolution to yield high-fidelity loads on very small components, such as the various protuberances and appendages 
on the Ares launeh vehiele. The Ares and Orion flows are inherently complex and can manifest nonlinear effects due 
to shocks, separated flows, strong jet interaetions, fluid/strueture interactions, and blunt-base wakes. Most of these 
flow features can only be efficiently grid resolved using solution-adaptive grid capability. The following sub- 
sections present some of the key enhancements to the grid generation capability that enabled the rapid generation of 
high-quality tetrahedral grids and flow solutions. 

A. NURBS Patches 

NURBS patches are a new versatile patch type incorporated in VGRID based on geometry routines from 
GridTool. NURBS patches can be defined with highly irregular boundary shapes and can be composed of an 
arbitrary number of sides. VGRID NURBS patches require an underlying NURBS surface definition. NURBS 
patches lead to a significant reduction in time spent for generating patches, because a NURBS surface definition can 
almost always be directly turned into one patch. The older parametric patches usually required extensive work to 
split a surface definition into 'good' patch shapes. Another benefit of NURBS patches is that grid points generated on 
NURBS patches are guaranteed to be on the original surface definition. Consequently, there is no need for time- 
consuming projection of the surface grid onto the underlying surface definition. Also, NURBS patches can 
accommodate multiple-loop patches, where a patch has one or more holes in the middle. This allows for an easy 
insertion, deletion or modification of small geometric features that may be embedded in a larger feature, such as a 
protuberance on a fuselage. These features can be patched and sourced in a separate file, and subsequently placed 
anywhere on the original configuration. The features may even be transformed and translated to different positions 
and orientations, and inserted in multiple locations. 

B. Volume and Surface Sources 

A new approach for the distribution of grid points on the surface and in the volume has been developed and 
incorporated in VGRID^^. In addition to the point and line sources, the new approach uses surface and volume 
sources for automatic curvature-based grid sizing and convenient point distribution in the volume. A new 
exponential growth function produces smoother and more efficient grids and provides superior control over the 
distribution of grid points in the field. Anisotropic grid stretching is still available for all types of sources for grid 
economy. The new approach makes use of an efficient bounding-box auxiliary medium for storing grid parameters 
defined by surface sources instead of three-dimensional background grids. These enhancements make the new 
approach less memory intensive and more computationally efficient. The latest version of VGRID has been used to 
generate a grid consisting of one billion tetrahedral cells around a generic wing-body transport configuration on an 
intel-Mac desktop. 

Figure 3 demonstrates an example of a grid around Apollo Launch Abort System using the enhanced VGRID. 
This configuration features several geometric complexities that present formidable challenges to grid generation 
such as truss beams in very close proximity with tight angles at end junctures, and grid extending through the abort 
motor nozzles and into the plenum. The grid in figure 3 was created from only seven volume sources around the 
vehicle and two volume sources in base wake region. The resulting grid of 48,355,970 tetrahedral cells and 
8,216,471 nodes was generated in 88 minutes of Central Processing Unit (CPU) time on a MacBook Pro laptop with 
a 2GHz Intel Core Duo processor and 2GB of RAM. The grid in the complex truss region is controlled by a single 
conical volume source that results in a high-quality thin-layered tetrahedral grid suitable for viscous flow 
simulations. 

GridTool (used to prepare the input file for VGRID) has also been upgraded with several usability 
enhancements and includes support for the latest VGRID (version 4.1) features. In addition, it has been ported to 
several compute platforms such as Mac and Linux. 
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C. Improved Gridding Practices 

Grid spacing around cylindrical features of a slender eonfiguration ean be preseribed with either a ring of line 
sourees or a single eylindrieal volume souree. Although a ring of line sourees ean add flexibility by allowing for 
loealized ehanges in grid spaeing, they also introduee undulations in the mesh as these sourees deviate from an exaet 
eirele. For the small eylindrieal features on the Constellation vehieles, sueh as extruding rings around the launeh 
vehiele or steps around the launeh abort vehiele tower, it was determined that eylindrieal volume sourees produee 
markedly better quality grids. The numerieal resolution of sueh features requires a high longitudinal eoneentration of 
eells with relaxed eoneentration in the eireumferential direetion. It was found that a eireumferential stretehing of 
10:1 was generally aeeep table, although the stretehing was lowered to 5:1 in some regions. 

Nozzle jet flow regions were resolved with approximately 50 eells aeross the plume. Best results were aehieved 
with uniform mesh size within the nozzle, whieh were generated with no stretehing using a truneated eonieal volume 
souree (a variation of the eylindrieal souree). Outside of the nozzle, the plumes were resolved by solution adaptive 
gridding, whieh offers more optimal grid point distribution. 

D. Solution Adaptive Gridding 

Solution adaptive gridding was aeeomplished by the CraftTeeh CRISP CFD® eode^^. CRISP CFD® is a stand- 
alone tool for mesh modifieation and quality improvement of three-dimensional mixed-element unstruetured 
meshes. Meshes eomprised of tetrahedral, prismatie, and hexahedral regions may be readily modified to generate 
more aeeurate flow solutions through loeal refinement and eoarsening. The adaption is driven by either an estimate 
of the solution error, or gradients of flow features. For tetrahedral eells as used in the present work, loeal refinement 
is aehieved by means of a eonstrained Delaunay refinement algorithm eombined with a eireumeenter point 
plaeement strategy. Any ineonsisteney between the eireumradius of a tetrahedron and some desired point spaeing 
triggers the point insertion proeedure. This iterative eell refinement is repeated until the eell eireumradii are 
eonsistent with the preseribed point spaeing. Coarsening of the tetrahedral region is also permitted through an edge 
eollapse proeedure. In regions where the grid is distorted or where solution errors are negligible, edges may be 
seleeted for removal. All eells ineident to the deleted edge are removed from the mesh, the adjaeent eells are 
redefined, and the two nodes of the edge are eollapsed to a single vertex. Capability exists to modify the viseous 
layers of tetrahedral grids and underlying surfaee mesh. However, our initial experienee has led us to only adapt in 
the inviseid region of the grid at this time. 

IV. Enhancements to Flow Solver 

The transition of USM3D from a researeh-oriented eode to a produetion-oriented eapability within a short time 
frame presented a ehallenge. The applications for the Constellation program required high throughput on very large 
grids. The immediate need for USM3D was to overeome several barriers to large-seale applieations sueh as memory 
bottleneeks, slow run times, and erippling robustness issues for high-Maeh number flows. Furthermore, 
improvements were required in the jet modeling eapability within USM3D to address inadequate jet boundary 
eonditions, diffieulties with starting jet flows, and defieieneies of turbulenee models for strong jet/aero -interaetion 
flows. These defieieneies prevented any meaningful use of the tool within this projeet environment. 

The following enhaneements were made to the flow solver to resolve these shorteomings. While largely driven 
by the urgent needs of the Constellation program, these enhaneements will elearly benefit the broader aeronauties 
eommunity. 

A. Overcoming Barriers 

7. Memory management for large-scale applications 

USM3D offers a large degree of flexibility in seleeting the number of proeessors for running a solution from a 
global grid and restart file. Earlier versions of the flow solver made extensive use of the eolleetive eommunieation 
paradigm and supported FORTRAN unformatted grid files with large reeords. However, these features earned a 
large memory footprint within the eode that impeded sealability of memory. The original implementation of 
parallel-proeessing inside USM3D was ideally suited for small to medium size grids (less than 20 million eells) 
running on shared memory maehines, sueh as NASA Advaneed Supereomputing (NAS) Columbia Altix. 
Unfortunately, this implementation was very ineffieient on the NAS Columbia maehine for grids over 50 million 
eells, due to exeessive swapping of memory buffers. 

The Constellation program presented a growing need to resolve intrieate geometrie details on a large vehiele as 
seen in figure 4. Consequently, grids grew from approximately 20 million tetrahedral eells to over 120 million 
tetrahedral eells. Attempts to run USM3D on sueh large grids at NAS were further hampered in 2008 by the advent 
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of distributed-memory maehines with smaller memory per eore (e.g., the Pleiades supereomputer with Xeon E5472 
Harpertown dual quad-core processors). To address this difficulty, several improvements were made to USM3D to 
reduce its runtime memory. For example, the enhanced flow solver could handle binary grid files and allowed for 
external pre-processing of grids for a parallel run. Additionally, extensive modifications were made to the flow 
solver data structures and inter-processor communication strategies to restore memory scalability. 

As a demonstration of this enhanced capability, a USM3D Navier-Stokes solution was generated on a civil 
transport configuration for the 4^^ AIAA Drag Prediction Workshop (DPW IV) using a grid of 212 million 
tetrahedral cells and the Spalart Allmaras turbulence model. This solution required a wall time of 4.62 seconds per 
iteration, using 700 grid partitions and 712 1-GB cores of NAS Pleiades Harpertown nodes. To date, the largest case 
to be successfully computed on the NAS Pleiades machine with the enhanced USM3D code is on a 623 million-cell 
mesh around the DPW IV geometry. 

2. Runtime Speedup 

The increased grid size also resulted initially in a significant increase in solution turnaround time, attributable to 
higher CPU and inter-processor communication times. In order to overcome this difficulty, USM3D’s internal grid 
pre-processing strategy was revised for a multi-processor run. Additionally, cache characteristics of the flow solver 
were improved. As a result of these enhancements, the speed up of the latest USM3D is five to six times the 2006 
version of USM3D. 

The upgraded USM3D has greatly increased its efficiency and productivity for the Constellation program. For 
example on an Ares I A106 task, the code was used to compute 800 high-quality Navier-Stokes solutions using a 90 
million-cell grid. These solutions were generated at a rate of 20 solutions per day using 6,000 CPUs on the NAS 
Pleiades cluster. 

3. Robustness for High-Mach Number Flows 

The Apollo Faunch Abort System was selected in the precursor studies for developing simulation guidelines for 
the Orion project FAS database. The intricate details in the Apollo FAS including truss geometry (see fig. 3) 
coupled with high Mach flows adversely impacted robustness of USM3D solutions and brought the FAS simulation 
guideline study to a halt. 

To address these difficulties, three new inviscid flux schemes, namely, HFFC, AUSM+, and FDFSS were 
implemented in the flow solver. Additionally, several enhancements were made to the implicit solution procedure. 
These improvements made it possible to readily obtain converged Apollo Faunch Abort System solutions up to 
Mach 6 using SA and SST turbulence models, as illustrated in figure 5. 

B. Modeling of Jet Flows 

The Ares series of launch vehicles utilize large numbers of control jets during ascent and stage separations to 
maintain stable flight. Similarly the Orion FAS utilizes powerful Abort Motor (AM) jets and smaller Attitude 
Control Motor (ACM) jets during a launch abort that dominate the aerodynamic characteristics of the vehicle. All of 
these jets operate at very high velocities and temperatures in close proximity to aerodynamic surfaces over a wide 
range of flight Mach numbers, from subsonic to supersonic. Within such flows, large gradients in pressure and 
temperature have a profound effect on the development of turbulence for transonic and supersonic jets. For high- 
speed, high-temperature jet flows, standard turbulence models lack the ability to predict the experimentally observed 
mixing rate of a shear layer, which affects plume spreading and results in large predictive errors for the 
aerodynamics of jet/surface interactions. Prior to the Constellation program, USM3D had limited capability for 
prescribing multiple jets and was deficient in the turbulence models for jet flows. Hence, the following 
improvements were made. 

7. Generalized supersonic boundary conditions 

The baseline USM3D flow solver allowed modeling of up to four nozzles using a cumbersome input procedure. 
In order to increase the productivity of USM3D, a generalized supersonic inflow boundary condition was 
implemented that not only allowed for an unlimited number of nozzle jets but also obviated the need for painstaking 
specification of nozzle orientations. These boundary conditions are prescribed on selected patch surfaces for each jet 
using a simple utility code that generates an additional input file for the flow solver. The approach to defining the jet 
boundary conditions for the applications to Ares launch vehicles and Orion FAS configurations is described in 
Appendix A. 

2. Initialization of jet flows 

Jet plume computations for Ares and Orion flight vehicles involved very high nozzle pressure ratios that 
presented severe startup problems with USM3D. In the earlier USM3D code, cells could only be initialized to 
freestream conditions at startup. When a jet boundary surface supporting a very high density, pressure, and 
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temperature falls adjaeent to volume cells initialized to freestream conditions, highly unstable and volatile non- 
physical numerical instabilities are produced. This was resolved by initializing a subset of cells inside the nozzle or 
plenum near the jet boundary surface with the approximate jet conditions. A simple scheme was contrived for 
initializing volume cells using the existing minimum distance function devised for the turbulence models. During 
preprocessing of the grid, cells within a desired distance from selected surface patches inside the nozzle or plenum 
are flagged for initialization. 

The robustness of this technique is demonstrated in figure 6 on the Apollo Launch Abort System. With a 
challenging nozzle pressure ratio of 4,200 in a freestream Mach number of 4, the solution readily started and 
subsequently converged into a strong jet/surface interaction flow. 

3. Extended 2-equation turbulence models 

The need for better physical modeling of jet flows in the USM3D turbulence models became apparent early 
during the Orion and Ares campaigns. Several approaches to address this deficiency have been advanced in the 
literature to modify the turbulence transport equations in a variety of ways. Some of the most proven and simple 
modifications to two-equation turbulence models have been selected and implemented in the USM3D flow solver 
Table 1 shows the range of turbulence models in USM3D. As denoted by cross-hatching in table 1, the 
modifications include the addition of correction terms for compressibility (Comp) and pressure dilatation terms (PD) 
in the turbulence transport equations for high-speed flows, and the addition of a simple modification to the 
Boussinesq’s closure model coefficient for high- temperature (Temp) jets. Two compressibility correction 
formulations are available, one from Sarkar and another from Wilcox. An extensive guideline study, yet to be 
published, was performed on the Orion LAS with Abort Motors and Attitude Control Motors that led to the selection 
of the SST with Comp+PD corrections as the preferred model for the jet- interaction computations. 


Table 1. Available turbulence models in USM3D ver. 20100611. 
(New enhancements^^ denoted by cross-hatching) 



V. Enhancements for Large-Scale Application 

A. Script-based run processes 

The effective use of TetrUSS in a production role required more than just improving VGRID and USM3D 
capabilities. Such tasks necessitate the rapid generation of dozens of grids and hundreds of flow solutions, which 
inherently present challenges in case management, data management, data analysis and post-processing, and quality 
control. These are process-related functions that were addressed using various scripting procedures. Scripts were 
developed to execute and manage processes for Ares and Orion, such as directory management, job submittals, best 
practices, post processing, data extraction, data plotting, and archiving. Examples will be included in the section on 
sample applications. 

B. Data analysis 

The Ares and Orion projects had some specific data analysis requirements that necessitated new processes and 
some modest development of capability. Process scripts were setup to 1) extract sectional load distributions, 2) 
subtract surface pressure values to create ACp contours plots, 3) a cross-sectional cutter utility for creating Cp line 
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plots, and 4) enhanced component force and moment integration capability to quantify loads on small features of the 
Ares rocket. 

1. Sectional Load Distributions 

The Ares structural and aeroelastic analyses require the aerodynamic loads to be defined at the launch vehicle 
centerline. A process was developed to transfer aerodynamic loads from a general three-dimensional unstructured 
surface mesh to the launch vehicle centerline. The process is based on the discrete data transfer technique developed 
by Samareh^^. This technique guarantees conservation of aerodynamic forces and moments to within machine 
accuracy. The algorithm is described in appendix B. 

2. ACp surface contour plots 

Jet interference effects were best discerned using ACp surface contour plots. This involved extracting the 
difference in contours of surface pressure coefficient, C^, between jets on and jets off conditions. A simple ACp 
surface contour process was developed using the TecPlot® analysis tool. The result is a color contour map displaying 
aerodynamic influence from the jets. Examples of this capability will be included with the sample applications in the 
next Section. 

3. Cross-sectional cutter for Cp line plots 

Capabilities within TecPlot® were exploited to automatically extract surface line plots of Cp. Process-control 
scripts were developed to create multiple comparison plots with experimental data. Examples of this capability will 
be demonstrated on the Ares and Orion vehicles in the next Section. 

4. Enhanced component force and moment integration 

An effort was initiated during the course of Ares I AD AC-2 A (AlOl) computational analysis to establish a 
procedure to isolate the contribution from each protuberance toward vehicle’s total integrated aerodynamic forces 
and moments. A sample result from this analysis is presented in figure 7 that displays axial force contribution from 
either an individual or a group of protuberances as a fraction of vehicle’s total axial force at various Mach numbers. 
The computed results in figure 7 readily indicate that among all protuberances, the Booster Deceleration Motors 
(BDMs), stiffeners/kick rings, and EAS nozzles principally contribute toward the total axial force. It can also be 
seen that the axial force contribution from BDMs peaks at Mach 0.9 and decreases thereafter with increasing speed. 
Similar analysis has also been performed for all other computed coefficients. Such a data analysis has aided the 
design of vehicle protuberances from an aerodynamic perspective. 

VI. Sample Applications 

Since 2006, an extensive amount of vetted computational aerodynamic and loads data have been delivered to the 
respective projects for the Ares I, I-X, and V launch vehicles and the Orion CEV and EAS configurations. Over 
7,000 USM3D solutions have been generated on the Ares vehicles, and well over 1,000 solutions produced on the 
Orion configurations. The following describes some typical applications of the TetrUSS toolset to the Ares and 
Orion vehicles. 

A. Ares I, I-X, and V Launch Vehicles 

7. Ares I Ascent Aero 

The Ares I Design Analysis Cycle (ADAC) outer mold line (OML) definition has evolved since 2006. During 
this evolution, aerodynamic databases have been developed for at least five major Ares I configurations with wind 
tunnel testing and CFD. CFD has played an important role in providing data for vehicle aerodynamics during ascent 
and first-stage-separation phases of a flight. CFD has also provided information on aeroelastic characteristics during 
the ascent phase. CFD solutions have been obtained using three flow solvers, namely, USM3D, OVERFLOW, and 
FUN3D. In the first two years of this project, extensive CFD computations were made to establish confidence in the 
ability of CFD to accurately simulate configurations similar to Ares I. This was accomplished via code-to-code and 
code-to-wind-tunnel comparisons. USM3D was selected as the primary flow solver for the aerodynamic assessment 
of Ares I during the ascent phase of flight. OVERFLOW was used to analyze aerodynamics of first-stage separation, 
whereas FUN3D was used to analyze aeroelastic characteristics of the vehicle. USM3D results have played a multi- 
faceted role in the Ares I ascent aerodynamic studies as described below. 

Due to inherent limitation of wind tunnel measurements associated with a slender geometry, comprehensive 
measurements of surface pressures is very difficult. Therefore, USM3D CFD results were anchored at selected 
locations and conditions with experiment and then used to provide more detailed vehicle surface pressures and 
sectional line loads. These data have been used by various groups for venting and debris analysis and load analysis 
for the Orion capsule and individual launch vehicle protuberances. Since early 2007, majority of USM3D results 
were completed prior to any wind tunnel testing. These predicted results aided wind tunnel model design, 
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protuberance fabrication, test matrix development, and instrumentation in the region of major flow variations. The 
USM3D predictions were validated against subsequently available wind tunnel test measurements^ As shown in 
the following, and in reference 24, good correlation between USM3D predictions and wind tunnel data further added 
confidence in the validity of the CFD results. Consequently, USM3D results were utilized in a number of trade 
studies to assess aerodynamic effects due to geometrical modifications and repositioning of the existing 
protuberances and addition of new ones to improve vehicle performance. 

The Guidance, Navigation, and Control (GN&C) aerodynamic database was primarily developed from the wind 
tunnel data on the 1 -percent model due to its ability to generate a large quantity of data in a timely manner. The 
USM3D force and moment results were utilized to provide flight Reynolds number increments to the wind-tunnel 
derived database as described in reference 25. 

USM3D was also used to generate the aerodynamic database for jet interactions from the Roll Control System 
(RoCS) on the Ares I full protuberance configuration. All USM3D computations were based upon ideal gas 
modeling of the flow. The baseline computational grid for this study was constructed following the best practices for 
the Ares I ascent power-off condition. Approximately 2 million grid cells surrounding each of the 12 RoCS thrusters 
were added to provide adequate grid resolution to support plume development and interactions with the free-stream 
flow. Reference 26 provides more details about this study. 

During the entire Ares I project, USM3D solutions for the ascent phase of the vehicle were based on tetrahedral 
grids that ranged from 1 0 million cells for the clean configuration at wind tunnel conditions to 90 million cells for 
full protuberance RoCS configuration at flight conditions. A converged solution typically needed about 20,000 
iterations and 5,000 CPU hours on NAS Pleiades Harpertown cluster. Nearly 1,000 USM3D solutions were 
generated within two months for the latest Ares I configuration (AD AC-3 A 106) to satisfy all project requirements 
for CFD data. Extensive use of scripts was made to set up (1) input stream for various freestream Mach numbers (2) 
grid rolling necessary for generating solutions at various roll angles (3) run directories for various solutions (4) user 
alerts for solution status, and (5) data extraction and graphical reduction of results. Use of scripts resulted in hands- 
off production of data in support of project requirements. Sample validation results from these studies are presented 
below. 

Figure 8 demonstrates convergence history for the Ares I A106 USM3D solution obtained at Mach 1.6, an angle 
of attack of 8^", a roll angle of 0^", and at wind tunnel Reynolds number. All the six components of forces and 
moments around the vehicle geometry were converged to within a certain percentage of their average values over 
the preceding 2,000 iterations. A value of 0.1% was used for the axial as well as normal forces and pitching moment 
whereas a value of 0.5% was used for the side force and rolling as well as yawing moments. As seen from the figure, 
lateral forces and moments took longer to achieve prescribed convergence criterion. The results shown in the figure 
8 have been automatically generated using scripts. 

The Ares I ADAC-2B (A103) configuration has been tested at three wind tunnel facilities^^'^^ covering low 
subsonic to high supersonic flow conditions and low to high Reynolds number. A sample comparison of a USM3D 
solution with wind tunnel measurements for this configuration at Mach 1.6, an angle of attack of T, and low 
Reynolds number wind tunnel conditions is presented in figure 9. Comparison of computed and measured surface 
pressures along the top centerline and around one of the largest protuberances (LH2) on the vehicle is presented in 
the figure. USM3D validation results for the same configuration using transonic high Reynolds number wind tunnel 
test^^ data can be found in reference 24. Assessment of USM3D results at flight Reynolds number is presented later. 

Comparison of the sectional normal load distributions over the Ares I A103 clean configuration as obtained from 
USM3D and wind tunnel data at various Mach numbers is presented in figure 10. As is evident from the figure, 
USM3D predicted the sensitivity of sectional loads to the changes in free stream Mach number. Sectional load 
extraction utility described earlier was used to derive sectional load variations along the stream wise direction. 

2. Ares I-X Flight Test Vehicle Configuration 

The Ares I-X is the first Flight Test Vehicle (FTV) for Ares I. For the Ares I-X FTV, USM3D flow solver was 
used to conduct all pre-flight simulations that were incorporated into the forces and moment databases. The merged 
liftoff, ascent and high angle-of-attack aero and line loads databases provide the basis for the GN&C algorithms and 
were used to calculate the trajectory and the footprint for the recovery of the first stage as well as for the loss of 
control or abort cases. 

The baseline configuration for Ares I-X FTV was Ares I AD AC-2 A (AIOI) full protuberance configuration for 
which a wind tunnel derived database existed. USM3D solutions were used to account for the changes in 
configuration between AlOl and Ares I-X FTV and to generate an Ares I-X FTV database at wind tunnel 
conditions. Additional USM3D solutions were generated at Ares I-X FTV flight conditions that provided Reynolds 
number increments between wind tunnel and flight. The latter increments were included in the dispersion bounds 
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that were provided to the GN&C team. The RoCS Jet Interaetion (JI) flight database was developed from hundreds 
of USM3D solutions. 

The generation of Ares I-X FTV aerodynamie database, and eomparison of pre-flight eomputations and flight 
data are detailed in referenee 27. A sample result from this effort is presented in figure 11 that shows longitudinal 
surfaee pressure distributions near the front of the vehiele at four azimuth loeations. The USM3D predietions were 
obtained using two different turbulenee models, namely, SA and SST. The eomputed pressures from two models are 
nearly identieal. The results show a very good eorrelation between USM3D predietions and flight measurements. 

Over the entire Ares I-X projeet, more than 5,000 USM3D solutions have been generated based on grids that 
eonsisted of 70-140 million tetrahedral eells. All of these solutions have been generated on NAS Pleiades and Altix 
supereomputers. A typieal Ares I-X solution was started and run with spatially first-order differeneing for a 
preseribed number of iterations, then switehed to globally seeond-order differeneing for the remainder of the 
eomputation. With most numerieal flow solvers, strong shoeks and sharp eorners ean often induee loeally non- 
physieal phenomenon sueh as spurious entropy, negative pressures, or exeessive pressure levels. For the Ares I-X 
eomputations, eells with pressure dropping below a minimum threshold of 0.0014 of freestream pressure or 
exeeeding a maximum threshold of 20-pereent above freestream total pressure were loeally switehed to first-order 
differeneing to eushion any potential instability. 

3 . Ares V Cargo Launch Vehicle 

The Ares V launeh system is in an early stage of planning and design with attendant revisions to vehiele 
eonfiguration. Due to limited availability of wind tunnel data, CFD solvers are being utilized to generate initial 
aerodynamie databases for the Ares V projeet. In this effort, OVERFLOW is the primary flow solver whereas 
USM3D is lending a supporting role providing approximately 40% of data. Referenee 28 provides details of 
database development and eomparisons of CFD and wind tunnel results. 

USM3D is also being used to establish best praetiees for investigating Plume Indueed Flow Separation (PIFS) 
eharaeteristies along the future Ares V heavy lift vehiele. In this preeursor study^^, USM3D is used to eompute PIFS 
along the Saturn V heavy lift launeh vehiele with F-1 roekets firing. Flight test PIFS data at four supersonie Maeh 
numbers are available for this eonfiguration. The USM3D solutions have been eomputed at the eorresponding 
supersonie freestream speeds, with angle of attaek of 0^", sideslip of 0^", and flight Reynolds number. In this study, 
sensitivity of the eomputed results to grid resolution, turbulenee models, and engine exit boundary eondition 
speeifieation approaehes is evaluated. As an illustration, the Maeh eontours along the Y=0 plane interseeting the 
Saturn V vehiele is shown in figure 12a at the Maeh = 6.5 eondition^^. PIFS ean be seen up the side of the vehiele 
from the nozzles by the blue eontours of low Maeh number flow. The Maeh eontours along the Y=0 plane 
interseeting the Ares V vehiele are shown in figure 12b for the Maeh 6.5 eondition. 

B. Orion CEV Launch Abort System 

The aerodynamie database development of the Orion Launeh Abort System (LAS) is led out of NASA Johnson 
Spaeeflight Center (JSC) under the CEV Aeroseienees Projeet (CAP) that provides multi-eenter eoordination for 
wind tunnel testing (WTT) and eomputational support. TetrUSS/USM3D is eontributing to CAP by investigating 
various sensitivities, sueh as WTT sting interferenee effeets, WTT-to-Flight Reynolds number sealing with Abort 
Motors (AM), AM thrust offset, and Abort-Motor/Attitude-Control-Motor (AM/ACM) jet interaetions. Eaeh of 
these tasks required a delivery of hundreds of flow solutions with stringent data requirements and very tight 
sehedules. Henee, the TetrUSS enhaneements diseussed earlier were eritieal for delivering a quality produet. The 
following will deseribe a typieal task within the CAP environment. 

A study was eondueted to quantify WTT-to-Flight Reynolds number sealing for the LAS 60- AA WTT 
eonfiguration with powered AM jets (fig. 13). This task required the delivery of 112 USM3D simulations of the 
Orion LAS with AM jets firing, in eonfigurations representing both free-flight and the 60- AA WTT. Data eovered a 
Maeh range from 0.7 to 2.5 and angles of attaek from 0^" to -Ih'’ in 4"^ inerements to mateh the available wind tunnel 
data. The objeetive was to produee data for evaluating USM3D flow simulations of WTT AM jet inerements and 
absolute magnitudes through eomparison with wind tunnel data and similar eomputations performed using Cart3D 
and OVERFLOW. In addition, the analysis of WTT and Flight eonditions provided WTT-to-Flight inerements that 
were used for database eorreetions and uneertainty evaluations. As in all of the tasks, the eomputational data had to 
be post-proeessed and delivered under speeifie requirements for arehiving. This was aeeomplished using numerous 
seripting proeesses. 

Four baseline grids were generated for this task. Typieal surfaee triangulations are shown in figure 14. Two grids 
were required for the vehiele with WTT AM nozzle inner mold lines (IML), and two for the IML of the Flight AM 
nozzle. Based on prior guideline studies, a finer baseline grid of approximately 25 million eells was required for the 
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transonic Mach conditions of 0.9 and 1.05, and a coarser grid of approximately 20 million eells was needed for all 
other Maeh numbers. A erosseut of the volume eells of the WTT eoarse grid is shown in figure 15. Note that the 
eells tend to be smaller elose to the geometry and eoarsen in the field away from the surfaee, whieh is eharaeteristie 
of the Advaneing Front Method. This grid does not provide adequate resolution to model the AM jet plumes and 
vehiele base flows. Henee, the CRISP CFD® grid adaption eode was used to refine the grid around errors of density, 
veloeity, and pressure as illustrated in figure 16 for Maeh 0.9 with the AM jet firing with a thrust eoeffieient, CT=4. 
Note the signifieant elustering of field eells in the AM jet plume and vehiele wake regions. For this ease, the 25 
million eell baseline grid was adapted resulting in an inerease to 34 million eells. 

Prior to the Orion LAS AM and ACM jet tasks, an extensive guideline study was performed to develop best 
praetiees and estimate solution uneertainty. A wide range of strategies for ereating the baseline grids, adapting the 
grid, and turbulenee models were explored. A representative analysis of the effeet of turbulenee model is presented 
in figure 17 for the LAS 60-AA eonfiguration with CT=4 and one adaptive grid strategy. This plot eompiles a 
weighted statistieal eomposite of how mueh eaeh solution deviated from the experimental data aeeording to: 

A = 2 \CxUSM3D,i /^X,exp,j “ l|/ S 

i=\ / i=\ 

where Cx~Cn} Cm s^nd Wi— 1 , 1 , 1 , 1 ,2 

The pitehing moment Cm was weighted twiee as heavily as the other eoeffieients beeause it earned more 
importanee and uneertainty in the aero database development. From figure 17, the SST Sarkar Comp + PD was 
seleeted as the turbulenee model of ehoiee having the least statistieal deviation from experiment. Note also the large 
deviations at Maeh 0.9. Wind tunnel Sehlieren imagery suggests that this deviation was most likely eaused by strong 
wall-interferenee effeets in the test data at transonie eonditions. 

The best praetiees developed from the guideline study were eaptured in seript-based run proeesses designed to 
exeeute and manage all runs. These seripts provided eomplete hands-off exeeution of preproeessing, job submittal, 
adapting the grid, terminating solution at eonvergenee, post-proeessing and plotting results, and arehiving the 
solutions. The following seript-eontrolled proeess was developed for LAS jet flow tasks: 

1. Run Ist-order USM3D flow solution on baseline grid for 3,000 iterations using SST Sarkar 

Comp + PD turbulenee model. 

2. Perform grid adaptation using CraftTeeh CRISP CFD® eode. 

a. Adapt to errors of density, veloeity, and pressure (p,u,v,w,p) 

b. Constrain minimum inviseid spaeing at 1.0 inehes (full seale) 
e. Interpolate baseline solution onto adapted grid 

3. Restart interpolated solution on adapted grid with Ist-order seheme for 1,500 iterations. 

4. Restart Ist-order solution with 2nd-order seheme. 

5. Automatieally stop 2nd-order solution when the following eonvergenee eriteria are met. 

a. Variation of Cl, Cd, Cm, Cn, and Ca less than 1-pereent over last 1,000 iterations 

b. Average of last 1,000 iterations within 1-pereent of average of preeeding 1,000 iterations 

e. Absolute differenee in the minimum and maximum values of the eoeffieients over last 1,000 
iterations less than 0.004 

d. Residual errors deerease by 3 orders of magnitude and remain eonstant or keep 
deereasing 

Extensive analysis and eonvergenee plots were generated during post-proeessing. Figure 18 shows a typieal 
erosseut of Maeh eontours with AM jets on and off, along with the eorresponding surfaee Cp eontours. Also shown 
is a ACp=CpjM_ON - Cp^AMjDFF surfaee eontour plot whieh illuminates the jet interaetions. Similarly, Cp line plots 
were ereated and eompared with WTT data as represented in figure 19. 

The adapted grids in this study ranged from 26- to 30-million eells, with the exeeption of the grids for transonie 
Maeh 0.9 and 1.05 that ranged from 32- to 35-million eells. The solutions were run on the NAS Columbia Altix 
eomputer typieally with 128 proeessors. Resouree usage ranged from 1,000- to 3,800-CPU hours per solution, with 
an average of 2,100-CPU hours per solution for the 112 eases. 
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VII. Concluding Remarks 

The NASA Constellation program is relying heavily on Computational Fluid Dynamie (CFD) tools that 
eomplement ground-based wind tunnel testing by aiding the design and aerodynamic database development for the 
Ares I Crew Launch Vehicle (CLV), Orion Crew Exploration Vehicle (CEV) and Eaunch Abort System (EAS). The 
TetrUSS system serves as the lead CFD tool within the Ares project, providing over 7,000 solutions to develop the 
computational aerodynamic database and loads for the vehicles up through supersonic speeds. The USM3D solver is 
also lending critical support to the aerodynamic development of the CEV and its EAS configuration with over 1 000 
solutions, but only after the implementation of considerable enhancements. 

The fast-paced application-driven environment has presented many challenges and opportunities to implement 
required capabilities within the TetrUSS system. Key enhancements added to the VGRID grid generator and its 
application practices enabled the routine generation of numerous viscous tetrahedral grids between 30 and 100+ 
million cells on the EAS and Ares flight vehicles. A commercial tool for autonomous adaption of the grid to the 
flow solution was seamlessly folded into processes for computing strong jet interaction flows on the LAS. Memory 
bottlenecks that limited USM3D application to grids in the 20 to 30 million-cell ranges were removed from the flow 
solver to extend its range to grids up to 623 million cells. Improvements were made to coding practices that led to 
code execution at least 5- to 6-times faster than a 2006 version. Robustness of computing the strong jet interactions 
of the LAS was enabled by the implementation of new flux schemes, better linearization, and new initialization 
functions. The accuracy of the jet interaction flows was improved by the implementation of compressibility and 
temperature corrections within several two-equation turbulence models. Scripting processes were developed to 
execute and manage thousands of flow solutions. Data analysis capability was developed for extracting section load 
distributions, ACp surface contour plots, and organizing extensive component loads. 

Over the past five years, this tool enhancement effort has largely been successful because of the close day-to-day 
collaboration between in-house code experts and the project application experts. The impact of the enhanced 
TetrUSS system on the Constellation program has been substantial. Although largely driven by the urgent needs of 
the Constellation program, these enhancements will clearly benefit the broader aeronautics community. 

Appendix A - Method for Prescribing Jet Boundary Conditions 

The methodology for prescribing the Ares and Orion LAS jet boundary conditions in ideal-gas flow solvers 
presented below was provided via. Excel spreadsheet by Mr. James Greathouse at NASA Johnson Spaceflight 
Center. The jet boundary conditions (BC) are defined on a plane within the nozzle slightly down stream of the throat 
as illustrated in the sketch below. A general assumption of one-dimensional isentropic flow is applied to determine 
the flow conditions at the nozzle exit plane, and then transfer them to the BC plane. 


BC plane, 




Exit 


Schematic of nozzle jet boundary condition prescription. 


The relations for one-dimensional isentropic flow through a nozzle are 
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Ptot 
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Ptot _ ^tot 


(A3) 

(A4) 


The following process is used to determine the jet boundary conditions on the BC plane. 
1 . Prescribe initial conditions and parameters: 

a. Freestream conditions 

M„,p„,T„, p„,q„, a^= Sref, Yair. and 

b. Jet plenum gas properties 

Yjet, Rjei , and Ttotjet 


C. Geometric areas for nozzle exit, throat, and BC plane 

Ae, At, and^sc 


2 . 


Expand plenum chamber conditions to nozzle exit plane using jet gas properties 

a. Using the prescribed value of yjeh adjust nozzle exit Mach number, Mejet, in equation A1 to match the 
choked area ratio with the geometric area ratio, 

(Ae^A )jQi—AQ/Af 


^ This defines Mejet 

b. Compute isentropic relations at nozzle exit by substituting yjet and Mejet^^to equations A2-A4 

(Pto/Pe)jet (Tto/TQ)jQf (^Pto/ Pe)jet 


3. Calculate nozzle exit plane conditions assuming jet gas properties 

a. Define a target nozzle thrust coefficient 

CT,target=targetjhrust/q^ref 

b. Adjust plenum chamber pressure, ptot in thrust equation 

Ct = VPtot liPtot !Pe) jet 'Ki^ + Y jeMl.jet) ' COS(25°)]/(g„5,^jr) 

until Ct matches Ct, target- The cos(25°) imposes the Cj of each offset abort nozzle (see fig. 16) in the 
body axial direction. 

^ This determines ptot^or the jet conditions 

c. Complete the nozzle exit conditions from known quantities 

Pejet ~P toJiP to/Pe) jet 
Te,jet~ Ttotl{Tto/Te)jet 
Pejet ~P e,jet ^(^jet Tejet) 

^e,jet~^e,jet ^^^^^(Yjet ^jetTeJet) 

4. Switch from jet gas to air and recalculate chamber conditions. The general assumption is to maintain the thrust. 

Since thrust is a function of exit-plane pressure and jet momentum, the values nozzle exit pressure, density, and 
velocity will remain unchanged. The primary impact of switching gas from jet to air is to recompute the nozzle 
exit temperature. This directly impacts the speed of sound and exit Mach number, which in turn alters the values 
of the isentropic relations in equations A2-A4 and the total conditions. 
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a. Adjust nozzle exit jet temperature and Maeh number for air 

Te,air' Pejet air Pejet) 

A/e, air ^ ejet ^ S(jrt( Yair ^air air) 

b. Assign remaining jet properties to nozzle exit as air 

Pe,air~ Pejet 
Pe,air ~ Pejet 
^e,air ~ ^ejet 

c. Perform sanity eheek of thrust by eomparing new Cr,air with Ct, target 

CT,air = \-Pe.air + Yair^lmj ' COS{25j]/{q^S,,f) 

d. Compute new isentropie relations at nozzle exit using equations A2-A4 with air properties Yatr and exit 
Maeh number, Me, air 

(P to/p e) air (Tto/Te)air (.Pto/ Pe) air 

e. Compute for total jet eonditions with assumed air 

Ptot,air~ Pe,air *(Plo/Pe) air 
Ttot,air~ Te,air ^(Tto/Te)air 
Ptot,air~ Pe,air^(Pto/Pe) air 


5. Caleulate A/A * at the BC plane just down stream of the throat 

a. Use equation A1 to eompute new (Ae/A'')air using Yair and Mg, air 

b. Compute new ehoked throat area A * is 

A Aq/ (Aq/A ) air 

6. Propagate the exit plane eonditions to the BC plane 

a. Using Yair, adjust Maeh number on the BC plane, MgQair, in equation A1 until its (A/A’')bc matehes Abc 

/A^ 

This defines MBc,air 

b. Substitute Yair and MBc,air into isentropie relations equations A2-A4 to eompute BC-plane ratios 

ip tot ^Pbc) air (Tfot I^bc) air i^Ptot ^PBc)air 

c. Compute the BC-plane properties 

PBC~Ptot,air^( Ptot ^PBc)air 
TBC^Tfot,air^ (Tto/TBc)air 
PbC~ Ptot,air(Pto/ PBc)air 
VBC=MBC^Sqrt( YairKirTBc) 

d. Compute USM3D non-dimensional BC values using freestream eonditions: 

Pusm3d~ PbC ^ P 00 
i^usm3d~ i^BC^ 

Pusm3d~PBC^( Yair Poo) 

7. Apply j et boundary eonditions 

a. Direetly apply pusm3d and Pusm3d aeross BC plane 
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b. Smoothly fan the velocity vectors from the center of the BC plane until parallel to the nozzle walls. 
Maintain nozzle axis component of constant across plane. 

c. Initialize the tetrahedral volume cells within the nozzle to approximate the conditions on the BC plane 
to facilitate startup. 


Appendix B - Launch Vehicle Line Load 

The Ares structural and aeroelastic analyses use line load coefficients that are defined as the axial derivatives of 
force coefficients (e.g., dC^/dx) on a set of user-defined points along the vehicle centerline. The launch vehicle 
surface is defined by a set of triangles, and CFD results are provided as pressure coefficients and skin friction 
coefficients at each node on the launch vehicle surface. The algorithm described here maintains conservation of 
forces and moments: the integration of line load will reproduce the CFD aerodynamic force and moment coefficients 
within machine accuracy. 

The line load coefficients are required at each user-defined point (z), which is surrounded by a control volume 
(Ax/) defined as: 




-,Ax^ =- 


3xi + X 7 X 

5i = = 


i-i+2xi + x,^i Xn_i + 3xn 

^ J A r — ^ 


where St represents the center of control volume surrounding node i (see fig. Bl). The integrated aerodynamic 
force and moment coefficients (AC) are calculated for each control volume (Ax^). This is a computationally intensive 
process that requires cutting surface mesh at discrete set of planes to calculate loads for triangles within the control 
volume and triangles partially cut by the control volume boundaries. The aerodynamic line load coefficients are then 
calculated as: 




CFD 



^1+5 


Figure Bl. User-defined point distribution and corresponding control 


Figure B2 shows a sample result for an Ares-I configuration. The blue line shows vehicle centerline with 
diamond symbols for the locations of user-defined points. Black line shows the outline of launch vehicle shape. The 
red line shows the variation of line load for normal force coefficient, Cn- 
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Figure B2. Sample line load coefficients 
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Launch Abort System 



a)- Major components. 


4 



b)- Nominal flight trajectory (from ref. 2). 


Figure 1. Major components and nominal flight trajectory of Ares-I and Orion flight 
vehicles. 
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Geometry Setup 
GridTool 



Grid Generation 
VGRID OpenGL 


TetrUSS is a loosely Integrated, flexible aero 

analysis and design system that can 

• provide initial solutions to complex problems 
within days 

• provide accurate representation of flow physics 

• be quickly adapted/extended to meet unforeseen 
challenges 

• have sufficient calibration and broad use for 
instiliing confidence 

• be easily applied by relatively non-expert users 



Flow Solver 
USM3D 





Visualization 
SimpleView 
(Commercial Packages) 


Tools & Utilities 


Figure 2. The Tetrahedral Unstructured Software System 
(TetrUSS). 
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2 Upper Stage Reaction Control Systems (2 US RCS) 

2 Upper Stage Roll Control Motors (2 RoCS) 

1 Umbilical Cover (1 Umbil) 

4 Upper Stage Ullage Motors (4 US Ullage) 

1 First Stage Tunnel (1 FSTUN) 

1 Upper Stage Tunnel (1 USTTUN) 

2 Booster Tumbling Motors (2 frustum BTM) 

2 Upper Stage Cameras (2 US CAM) 

1 Launch Abort System Tunnel (1 LAS tunnel) 

2 Upper Stage Ultra High Frequency Antennae (2 US UHF) 
8 Booster Deceleration Motors (8 BDMs) 

1 Upper Stage Liquid Hydrogen Feedline ((1 US LH2) 

4 Launch Abort System Nozzies (4 LAS Nozzles) 

1 Upper Stage Fill/Drain-Line (1 US F/D-line) 

4 Wedges 



Translucent Front View 



frustum region 


LH2 region 



CEV region 


Figure 4. Front view of ARES I ADAC-2B (A103) configuration with 
protuberances and a view of computational grid around selected region. 


Mach 4 Mach 5 Mach 6 





Figure 5. USM3D solutions for Apollo Launch Abort System using SA and SST 
turbulence models at free stream Mach 4, 5, and 6. Solutions obtained without 
plumes. 
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Figure 6. USM3D solution for Apollo Launch Abort System using SST turbulence model 
for free stream Mach 4, Nozzle pressure ratio 4,200. 
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Figure 7. Individual/group protuberance(s) contribution to the computed Cap for the 
ADAC-2A (AlOl) Cl configuration, a=7®, WT Reynolds number. 
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Figure 8. Typical solution convergence for the Ares I ADAC-3 (A106) configuration, 
M=1.6, a=8“, 0=0®, WT Re. 



Figure 9. Computed and measured surface Cp on the Ares I ADAC-2B (A103) 
configuration, M=1.6, a=7“, WT Re. 
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Figure 10. Computed and measured line loads over the frustum for the clean ADAC-2B 
(A103), a=7®, WT Re. 


Ares IX FTV Pressure Coefficient SA 

RoCS:cw FS Motor: on SST 

♦ Flight 






Figure 11. Comparison of USM3D surface pressures with the flight test data for 
the Ares I-X FUght Test Vehicle. Computations using SA and SST turbulence 
models (from ref. 27). 
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(a) (b) 

Figure 12. Symmetry plane Mach contours for two heavy lift launch vehicle 
configurations at M=6.5, a=0 , P=0 , flight Reynolds number, (a) Saturn V rocket (b) 
Ares V rocket (from ref. 29). 



Figure 13. The Project Orion LAS 60-AA wind tunnel model mounted in the 
NASA Ames Research Center 11-foot tunnel. 
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Figure 14. Coarse grid surface triangulation for LAS 60-AA configuration. 



Figure 15. Cross cut of baseline coarse grid (20 million cells) for the LAS 60-AA 
WTT configuration. 
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Figure 16. Solution adapted grid from guideline study on LAS 60-AA for 
Moo=0.9, a=-8 , AM_CT=4, SST Sarkar Comp + PD turbulence model (adapted 
grid has 34 million cells, baseline grid has 25 million cells). 


sum of mean weighted dev with CT=4 & BL1A1 



Figure 17. Statistical assessment of USM3D turbulence model variations from 
experimental data for LAS 60-AA. 
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Mach Contours 


Surface Cp and ACp Contours 


AM jet on 



AM jet off 
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Cutting Plane through Jets ((j) = 45'^) 


Side View 


Figure 18. Typical Mach crosscut and surface Cp, ACp contour plots for the LAS 60- 
AA configuration. M«j=2.5, a=-8 , ReD=2.5xlO^ AM_CT=3.689. 
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Figure 19. Topical surface Cp line plots for the LAS 60-AA configuration. 
Moo=2.5, a=-8°, ReD=2.5xlO^ AM_CT=3.689. 
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